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ABSTRACT 

Recent models of black hole growth in a cosmological context have forwarded a paradigm in which 
the growth is self-regulated by feedback from the black hole itself. Here we use cosmological zoom 
simulations of galaxy formation down to z = 2 to show that such strong self-regulation is required in 
the popular spherical Bondi accretion mod el, but that a plausible alte rnative model in which black hole 
growth is limited by galaxy-scale torques (jHopkins fe Quataertll201ll ) does not require self-regulation. 
Instead, this torque-limited accretion model yields black holes and galaxies evolving on average along 
the observed scaling relations by relying only on a fixed, 5 % mass retention rate onto the black hole 
from the radius at which the accretion flow is fed. Feedback from the black hole may (and likely does) 
occur, but does not need to couple to galaxy-scale gas in order to regulate black hole growth. We 
show that this result is insensitive to variations in the initial black hole mass, stellar feedback, or other 
implementation details. The torque-limited model allows for high accretion rates at very early epochs 
(unlike the Bondi case), which if viable can help explain the rapid early growth of black holes, while 
by z ~ 2 it yields Eddington factors of ~ 1-10%. This model also yields a less direct correspondence 
between major merger events and rapid phases of black hole growth. Instead, growth is more closely 
tied to cosmological disk feeding, which may help explain observational studies showing that, at least 
at z > 1, active galaxies do not preferentially show merger signatures. 

Subject headings: Black Hole Physics — galaxies: active — galaxies: evolution — quasars: general 



1. INTRODUCTION 

There is by now ample evidence that supermassive 
black hol es reside at the center of most sizeable galax- 
ies (e.g.. iFerrarese fc Fordl 120051 ) . During the last two 
decades, observations of local galaxies have established 
a number of well-defined relationships between the mass 
of the central black hole and stellar properties of the 
host galaxy such as the mass of the centr a l bulge (Mrh~ 
Mh„w relation; iMagorrian et al.l 119981 : lHaring fc Rix 



2004 ). the concentration or Sersic index (|Graham et al 
20011) the bulge velocity dis p ersion (Mbh~Q " rela- 
tion; IFerrarese fc Merrittl l2000l iTremaine et al.l 120021: 



iGultekin et all 120091: iMcConnell et alJl2011| ).~the virial 
bulge mass dMarconi fc Huntl [20031) . the spheroid bind- 
ing energy (jAller fc Richstonell200'7l). an d even the mass 
of the dark matter halo (Ferraresc 2002). The tightness 
of these correlations suggests that there is a close con- 
nection between the formation and evolution of galax- 
ies and their central black holes. Moreover, the simi- 
larity between the global cosmic star formation history 
(iMadau et al.l I1996T) and the ev olution of quasar abun- 
dances (|Bovle fc TerlevichlH998l) . as well as the apparent 
connection of starburst events and active galactic nu- 
clei (AGN) in individual objects (Kauffm ann et al.ll2003l) 
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provide further support for a link between galaxy forma- 
tion and accretion onto a central black hole. Understand- 
ing this connection can have broad implications in galaxy 
evolution and the formation of massive black holes in the 
early universe. 

Analytic arguments suggest that the observed black 
hole-galaxy correlations might arise from the coupling of 
energy and/or momentum from radiation emitted by the 
black hole accretion process to the surrounding galaxy (a 
process known as black hole or AGN feedback) , leading to 
the idea of self-regulated growth and the cq-evolution of 
black holes and galaxies (iSilk fc R ees 1998; F abianlll99S 
Kind 120031; [Begelman fc Nathl 120051: iMurrav et alJ koOE 



McQuillin fe McLaughlin! 120121) . The energy produced 



by accretion onto the black hole can easily exceed the 
binding energy of its parent galaxy and, therefore, even 
if a small fraction of this energy is transferred to the sur- 
rounding gas, black hole feedback can ha ve a profoun d 
effect on the evolution of its host galaxy ()Fabian|[2012l) . 
This interpretation has been particularly attractive for 
galaxy formation models since both semianalytic mod- 
els and cosmological hydrodynamic simulations require 
the injection of additional energy to prevent the over- 
cooling of ga s and the subsequent runaway formation 
of stars (e.g.. iCroton et al.l 120061: [Somerville et~aH 120081: 



IGabor et all 1201 It iTevssier et all 1201J. AGNfeedback 
is often invoked as the source of this energy. 

Motivated by these ideas, an increasing num- 
ber of studies have now incorporated models for 
black hole growth and feedback into idealized galaxy 
scale simul ations as well as fu l l cosmological simula- 
tions (e.g., iSpringel et all 120051: iDi Matteo et all 120081: 



Booth fc Schavel (20091: iDebuhr et all 120111: iKim et al 



2011 



2012). 



IPower et all 1201 lHChoi et al.l 120121 IDubois et aT 
Since the numerical resolution required to self- 
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consistently follow gas inflows from ~ lOkpc host-galaxy 
scales down to the central black hole (~ 10~ 6 pc for 
the Schwarzschild radius of a 10 7 M© black hole) is out 
of reach for current computers, simulations must rely 
on "sub-grid" models to infer black hole accretion rates 
based on the physical properties of the surrounding gas 
at the scales resolved in the simulation. The majority 
of galaxy formation simulations to date have employed 
black hole accretion prescriptions based on the spher- 
ical Bondi-Ho yle-Littleton parameterization (hereafter 
"Bondi rate" ; iHovle fc Lvttletonl 119391 iBondi fe Hovlel 
11944 iBondil |1952() , with the exception of some recent 
models that incorporate the physics of angular mo 
ment um transport in non-spherically symm etric flows 
(e.g., IDebuhr et all I20TI IPower et al.l l20Tl . A vari- 
ety of feedback models have been implemented to date, 
including the injection of a fraction of the accretion 
luminosity as therm al energy into the g as surround- 
ing the black hole (ISpringel et all 2005; S iiacki et al.l 



20091 ) . the injecti on of momen- 



l2007t Uohansson et aD _ 
turn from optically thick radiation fields (jDebuhr et al.1 
120 111) , heating and ra diation pressure from X-ray gen- 
erate d by the AGN (lHambrick et al.l 1201 It iKim et al.l 
120111 ; IChoi et all I2012D . and the injection of mass 
and / or momentum from AGN driven winds or radio- 
iets (IChoi et al.ll2012t IDebuhr et all 120121 ; iDubois et afl 
12011) . 

The success of coupled black hole accretion and feed- 
back models in accounting for a number of observa- 
tions including black hole mass functions, Eddington 
ratio distributions, and bla c k hole-galaxy correl ations 
(e.g., iDi Matteo et al.l 120051 : IHopkins et ahl 120061 ) have 
contributed to establish a paradigm in which black hole 
growth is self-regulated by feedback from the black hole 
itself. There is, however, no direct observational evidence 
to date for such self-re gulated growth of black ho les due 
to their own feedback | Alexander fc Hickoxll20r2l ). AGN 
feedback acting in the kinetic or radio mode is commonly 
observed in the form of radio-emitting relativistic jets, 
and its effects on the hot gas in galaxies and clusters are 
relatively well understood. The radiative or quasar feed- 
back mode has been identified by the presence of highly 
blue-shifted absorption and emission lines and broad line 
wings, but the overall impact of qu asar-driven w inds on 
the host galaxy remains uncertain ([Fabian 20"l2l). Even 
if AGN feedback is effective in quenching star formation 
in galaxies and suppressing cooling flows in groups and 
clusters (as required by galaxy formation models), there 
is a priori no reason for AGN feedback to be the dom- 
inant physical process regulating black hole growth in 
galaxies. 

In the presence of a non-zero net angular momen- 
tum of the surrounding gas, black hole feeding likely oc- 
curs throu gh viscous transport in a Keplerian accretion 
disk (e.g., IShakura fc Sunyaevl IT973I : iBalbus fc Hawlevl 
119981 ). Therefore, it is plausible that most of the feed- 
back energy and momentum escapes along the polar di- 
rectio n without significan tly affecting further accretion 
(e.g., lOhsuga et al.l l2005h . In order to allow for signifi- 
cant black hole growth, the accretion disk must be con- 
tinuously replenish by gas inflows from larger scales. The 
angular momentum of the gas along with competition 
with star formation are known to represent significant 



barriers to the transport of gas from gal actic scales down 
to the black hole accretion disk (e.g., iThompson et al.l 
120051; |jogeil2006l; IHopkins fc Quataertll2010D . Thus, the 
rate of angular momentum transport relative to the lo- 
cal star formation rate (SFR) may well be the dominant 
physical process governing black hole growth in g alaxies, 
even i f AGN feedback is acting at some level (e.g.. lEscalal 
I2006L 120071: IPower et al.ll2^lTf) . 

Numerical simulations have shown that large-scale 
tidal torques produced by galaxy interactions, galaxy 
mergers, and gravitational instabilities in self-gravitating 
disks can lead to efficient angular momentum trans- 
port and the rapid infl ow of gas into t he central kilo- 
parse c of galaxies (e.g.. lHernquistl lT989; Shlo sman et al.l 
Il989t iBarnes fc Hernquistl 119920 ! At smaller scales, 
large-scale torques become less efficient and additional 
mechanisms are required. Using multiple nested galaxy 
scale simulations o f prog ressively higher resolution, 
IHopkins fc Quataertl (|2010h showed that a series of grav- 
itational instabilities can generate net inflow rates of up 
to l-10M Q yr _1 down to sub-parsec scales (< O.lpc). 
They found that non-axisymmetric perturbations to the 
stellar gravitational potential produce orbit crossings and 
shocks in the gas that can be efficien t in removing angular 
mome ntum even at scales < 10 pc. IHopkins fc Quataertl 
(|2011[) derived analytic expressions for the loss of angular 
momentum in the presence of such shocks and presented 
an analytical "gravitational torque" model for the result- 
ing gas inflow rates. This analytic model was found to 
reproduce gas inflows rates at < 0.1 pc scales from simu- 
lations with significantly less scatter than the Bondi pa- 
rameterization for the total enclosed mass within a given 
radius. 

The goal of this paper is to investigate the growth 
of supermassive black holes due to gravitational-torque 
driven accretion in a cosmological context and deter- 
mine whether this growth mechanism can account for the 
observed black hole-galaxy correlations without the re- 
quirements of AGN feedback and self-regulated growth. 
Our approach consists on combining high resolution cos- 
mological hydrodynamic simulations together with ana- 
lytic models of black hole accretion in post-processing. 
This simplification allows us to follow the growth of 
galaxies and black holes from early epochs down to z = 2 
without making any prior assumptions about the effects 
of AGN feedback on galactic scales. By comparing pre- 
dictions from the Bondi and gravitational torque mod- 
els we evaluate the relative importance of feedback and 
angular momentum transport in regulating cosmological 
black hole growth. In addition, we discuss the implica- 
tions of different black hole accretion models on the link 
between major galaxy mergers and rapid phases of black 
hole growth, as well as the growth of black hole seeds 
in the early universe. The observed black hole-galaxy 
correlations are regarded here as a strong constraint for 
theories of black hole growth and galaxy evolution but 
not as a mandatory outcome of AGN feedback, as it is 
commonly assumed. Given that both star formation and 
black hole accretion histories peak at z ~ 2, our simu- 
lations target the critical epoch when such correlations 
were likely established. 

This paper is organized as follows. We describe our 
simulations and main analysis procedures in Section [2j 
In Section [3l we present and overview of our simulated 
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galaxies and describe their main physical properties. In 
Section|4l we evaluate the Bondi and gravitational torque 
models and show how the inferred black hole accretion 
rates relate to the physical properties and evolution of 
the host galaxies. In Section[SJ we use the observed black 
hole-galaxy correlations to put constraints on black hole 
accretion in cosmological time scales and evaluate the ef- 
fects of AGN feedback as required by different black hole 
models. In Sections [6] and [7l we show that the gravita- 
tional torque model naturally yields the observed black 
hole-galaxy correlations regardless of the initial masses 
of black holes and galaxies, and without the need for 
self-regulation by galaxy-scale AGN feedback. Finally, 
we summarize our results in Section [8l 

2. SIMULATIONS 

We use the ex tended version of the Gadget-2 code 
(ISpringell I2005D that has been fully described in 
lOppenheimer fc Pavel ([2008D . Gadget-2 combines a 
TreePM algorithm for computing gravitational forces 
with an entropy-conse rving formulation of smoothe d par- 
ticle hydrodynamics (jSpringel fc Hernquistl [20021 ). We 
include radiative cooling f rom primordial ga s assuming 
ionization equilibriu m as in iKatz et all (119961) and metal- 
line cooling based on iSutherland fc Dopital (|1993l) . Pho- 
toionization heating is modeled via a spatially uniform , 
optically thin UV background (lHaardt fc Madaul l20"0ll) 
starting at z = 9. 

Star formation and feedback are modeled in Gadget- 
2 through a sub-grid prescription. Gas particles that are 
sufficiently dense to become Jeans unstable are treated 
as a multi-phase fluid with col d, dense clumps embed- 
ded within a hot phase medi um (|McKee fc Ostrikerll977t 
iSpringel fc Hernquistl [20031 ) . Stars form from the cold 
phase and feedback from Type II supernovae causes evap- 
oration of cold clumps into the hot phase (which can con- 
dense back onto cold clumps) and metal enrich ment. The 
result ing SFRs are in accord with the observed lKennicuttl 
(|1998t) relation. We also account for energy and metal 
feedback from Type la supernovae and mass-loss and 
metal enrichment from asymptotic giant branch stars, 
tracking the production of four metal species (C, O, Si, 
and F e) separately as described in lOppenheimer fc Pavel 
(|2008l) . In addition, our simulations include a galactic 
outflow mechanism that imparts kinetic energy to gas 
particles. The outflow rates scale with galactic veloc- 
ity dispersion, a, and the ratio of material entering the 
wind versus forming stars (i.e., the mass loading fac- 
tor) scales as l/o", a s in the momentum-driven wind case 
([Murray et al.ll200"5l ). An on-the-fly galaxy finder is used 
to compute galaxy masses which ar e converted to a using 
standard relations (|Mo et al.lll998| ). 

Overall, this treatment of metal production and galac- 
tic winds has been successful in reproducing a broad 
range of observations including the chemical enrich- 
ment of th e IGM at z > 2 ([Qppenheimer fc Pavel 
120061 120081) the luminosity function of high-r edshift 
galaxies ([Pave et al.1 120061 : iFinlator etafl 120071) . and 
the galaxy mass- metallicity relation ([Finlator fc Pavel 
1200a iPave et all 120111 ). In addition, the high resolu- 
tion zoom simulations presented here have been suc- 
cessful in reproducing many of the star formation, mor- 
phological, and kinematic p roperties of z = 2 galaxies 
(|Angles- Alcazar et al.ll2013l ). 



We note that our simulations do not include prescrip- 
tions for black hole growth and the effects of AGN feed- 
back. Our goal is to use galaxy formation simulations to- 
gether with analytic black hole accretion models to study 
the physical processes governing cosmological black hole 
growth, based only on the observed black hole-galaxy 
correlations, and with no implicit assumptions about the 
amount and efficiency of feedback from the accretion pro- 
cess. 

2.1. Simulation runs 

In this work we use the "zoom-in" technique (e.g., 
INavarro fc Whit^[l99l to run cosmological simulations 
and resolve the physical conditions at the inner kilo- 
parsec of galaxies over cosmic time. We first ran an 
intermediate-resolution cosmological simulation with 2 x 
256 3 gas+dark matter particles in a [24 /i _1 Mpc] 3 comov- 
ing box. This simulation allows us to characterize the 
population from which galaxies are selected for resimula- 
tion. We chose 2 galaxies at z = 2 characterized by simi- 
lar masses but different morphologies, environments, and 
merger histories. We identified all particles within the 
virial radius of each galaxy and traced them back to their 
locations on the initial grid, where they mark the refine- 
ment region. Zoom initial conditions were then generated 
by populating the refinement region with a larger num- 
ber of lower mass particles (each particle is replace by 
8', where I is the zoom level), and adding the additional 
small scale fluctuations appropriate to the new Nyquist 
frequency. In order to reduce numerical artifacts due to 
the difference in particle masses and to make sure that 
the large scale gravitational torques acting on the target 
halos are accurately represented, the high resolution re- 
gion was significantly enlarged by an iterative "cleaning" 
procedure and surrounded by two nested, concentric lay- 
ers of progressively lower resolution. After running the 
simulations, we enlarged our galaxy sample by including 
6 additional galaxies that were found well within the high 
resolution regions, confirming that there is no contami- 
nation of low resolution particles within the virial radius 
of each galaxy. 

Our zoom simulations have an effective resolution 
equivalent to 2 x 1024 3 particles homogeneously dis- 
tributed in a [24 h~ 1 Mpc] 3 box, with (high resolution) 
gas particle mass m gas ~ 2.3 x 1O 5 M , dark matter par- 
ticle mass todm w 1.2 X 10 6 M Q and softening length 
e w 0.47/i _1 kpc comoving (or ~ 224 pc physical at 
2 = 2). All results presented in this paper (except where 
indicated) correspond to simulations including our pre- 
ferred model for galactic outflows, based on momentum- 
conserving scalings as explained above. In order to iden- 
tify the possible effects of winds in our results, we also 
ran all simulations without winds. Additionally, we ran 
simulations with 2 x 512 3 effective resolution (i.e., 2 times 
lower spatial resolution and 8 times lower mass resolu- 
tion) for all galaxies to test for numerical convergence. 

Throughout this paper we assume a ACPM cosmology 
with parameters f2 A = 0.72, Q M = 0.28, f2 b = 0.046, h = 
0.7, as = 0.82, and n = 0.96, cons istent with the five-ye ar 
WMAP concordance cosmology ([Komatsu et al.l[2009f l . 

2.2. Analysis of simulations 
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Figure 1. Snapshots of our sample of simulated galaxies at z = 2. For each galaxy, we show star formation rate surface density maps 
color-coded according to the projected density-weighted temperature (top panels) and the projected stellar surface density (bottom panels). 
The total star formation rates and stellar masses of galaxies are indicated in the top and bottom panels respectively. Face-on and edge-on 
views of galaxies correspond to the direction of the angular momentum at z = 2. The region shown is 30kpc (physical) across. 



We produced over 230 snapshot files at time inter- 
vals ranging from ~ 5 to 25 Myr during the simulation. 
We identify galaxies by means of the Spline Kernel In- 
terpolative Denmax algorith m (sKirfl) at eac h redshift 
snapshot independently (see lKeres et al.ll2005l for a de- 
scription). Our simulated galaxies are thus defined as 
bound groups of star-forming gas particles (i.e., gas par- 
ticles with densities above the threshold for star for- 
mation) and star particles. We initially associate each 
SKiD-identificd galaxy with a dark matter halo by us- 
ing a spherical overdensity algorithm, defining the virial 
radius as the radius enc losing a mean density given by 
iKitayama fc Sutol (|1996| ). Overlapping halos are then 
grouped together so that, by construction, every halo in 
our final catalog has a central galaxy (the most massive 
galaxy) and a number of satellite galaxies. Beginning 
at z = 2, we reconstruct the evolution of each individ- 
ual galaxy back in time by identifying each galaxy's most 
massive progenitor at all redshifts. All galaxies presented 
here are centrals. 

When calculating morphological and kinematic prop- 
erties of each galaxy as a function of time, we take the 
position of the most bound gas particle as the center of 
the galaxy, which corresponds to the position of the black 
hole. This results in a physically more meaningful and 
a significantly more stable definition compared to the 
galaxy center computed by skid, especially during close 
galaxy encounters and galaxy mergers. Since we are in- 
terested in the evolution of galaxies and black holes over 
cosmological time scales, we smooth the time evolution of 
all physical quantities and extract their main evolution- 
ary features by calculating averages over time intervals 
of 50 Myr (except otherwise noted) . Finally, we impose 
a resolution limit on our sample of galaxies by rcquir- 

1 http:/ /www-hpec. astro. washington.edu/tools/skid. html 



ing that they contain at least 100 gas particles and 100 
star particles within the inner kiloparscc. In this way, we 
ensure that the morphology and the physical conditions 
at the centers of galaxies can be accurately character- 
ized for the purposes of this work. Therefore, lower mass 
galaxies take part of the analyses and results only until 
they are resolved at later times. 

3. GALAXY PROPERTIES 

Our zoom simulations reproduce many characteristic 
features of observed z sa 2 star forming galaxies. The 
properties of these g alaxies have been analyzed i n detail 
in a separate paper (| Angles- Alcazar et al.ll2013T) . where 
we show that strong stellar winds are required in order 
to maintain high gas fractions, redistribute star-forming 
gas over larger scales, and increase the velocity dispersion 
of simulated galaxies, in better agreement with the large, 
extended, rotation-domin ated yet turbulent star-formin g 
disks of the SINS survey (jForster Schreiber et aT1l2009t) . 

Here, we provide a general description of all simulated 
galaxies and then we focus on the physical conditions at 
their centers, since they eventually determine the accre- 
tion rates onto the central black hole. 

3.1. Global properties 

Figure [T] gives a visual impression of our galaxy sam- 
ple by showing the distribution of the star-forming gas 
and stellar contents at z = 2. For most galaxies, the 
star-forming gas lies primarily in a thick, rotationally 
supported disk, with sizes ranging from ~ 5 to ~ 25kpc. 
A wide range of morphologies can be identified, including 
prominent spiral arms (galaxy g54) , a very compact gas 
distribution (galaxy g2403), a pair of interacting galaxies 
(g2438), and a turbulent, highly disturbed galaxy (g222). 
The projected stellar distributions also reveal a wide va- 
riety of galaxy sizes and morphologies. Most galaxies 
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Figure 2. From top to bottom: Stellar mass, gas fraction, and 
specific star formation rate as a function of time from z = 8 to z = 2 
for all simulated galaxies. Each color corresponds to a different 
galaxy. 

appear to have a large scale stellar disk surrounded by a 
significant spherical component. 

At z = 2, our simulated galaxies span about an order 
of magnitude in both stellar mass (from ~ 1.6 x 10 9 to 
- 3.3 x 10 10 M Q ) and SFR (from ~ 1 to ^Moyr" 1 ). 
Figure shows the evolution of individual galaxies from 
z = 8 to 2 in terms of their stellar mass, gas fraction, 
and specific SFR. Major mergers (identified as abrupt 
changes in the stellar mass of galaxies) can have a signif- 
icant impact on the dynamical properties of galaxies but 
represent only a small fraction of the total mass growth 
down to z = 2. Specific SFRs follow a common trend 
for all galaxies, decreasing from ~ lOGyr -1 at z = 8 
down to ~ lGyr -1 at z = 2, but their individual star 
formation histories are clearly different. Gas fractions 
evolve relatively slowly from nearly unity at z = 8 down 
to <~ 0.5 at z = 2, due to the early suppression of star 
formation by momentum-driven winds and the recycling 
of gas at later times. Overall, our simulated galaxies 
are characterized by a wide range of masses, morpholo- 



gies, and star formation his tories and are representativ e 
of "normal" z = 2 systems ([Angles- Alcazar et al.ll20 13h . 

3.2. The central kiloparsec 

In order to investigate the dominant processes govern- 
ing black hole growth, it is crucial to understand the 
physical conditions at the inner regions of galaxies, how 
they evolve in time, and what the main evolutionary trig- 
gers are. Here, we evaluate the masses and morphologies 
of the stellar and gas contents of galaxies as well as the 
specific SFRs within a radial distance Rq = 1 kpc from 
their centers. Large-scale gravitational torques produced 
by galaxy mergers and interactions are often invoked 
as triggers of star formation and AGN activity (e.g., 
ISpringel et al.lr2~005[) . The identification of galaxy merg- 
ers in cosmological simulations is, however, a non-trivial 
task in which the simple working definition of "galaxy" 
can result in a rather different timing of merger events or 
the temporal ident ification of close ga laxy encounters as 
a merging system (|Gabor et al.ll201lj) . For the purposes 
of this work, we simply associate changes in galaxy prop- 
erties with galaxy interactions and/or mergers by visual 
identification from gas and stellar surface density plots 
(similar to those in Figured]) at the appropriate snapshot 
times. 

Figure [3] characterizes the masses and morphologies of 
the stellar and gas contents of galaxies g54 and g222 
within the inner kiloparsec as a function of time. In the 
remainder of this section, we use these two galaxies to 
illustrate the analysis but we note that our conclusions 
apply to all simulated galaxies. Galaxy g54 undergoes a 
relatively quiescent evolution below z w 3, after a merger 
that results in a significant increase in the stellar mass 
within R a (~ 2000 Myr after the Big Bang). Galaxy 
g222 is located in a denser environment and it undergoes 
a more violent evolution, with mergers that result in the 
rapid increase of stellar and gas masses within Rq at 
times - 1500, 1900, 2700, and 3000 Myr. 

Here we parametrize the gas and stellar morpholo- 
gies within the inner kiloparsec of galaxies based on the 
amount of rotational support, K ro t; and the bulge mass 
fraction, /bulge- We calculate K ro t as the fraction of ki- 
netic energy in ordered rotation with respect to the tota l 
angular momentum within i?o fe.g.. [Sales et al.1 120101 ). 
Additionally, we estimate /bulge from a bulge-disk decom- 
position based on the full three-dimensional kinematic 
information of particles of either the gas or the stellar 
components. We calculate the azimuthal velocity (v^,) 
of each particle with respect to the direction of the to- 
tal angular momentum within Rq. In the presence of a 
spherical component supported by random motions, v<j> 
is expected to have a symmetric distribution with 50 % of 
the particles having velocities < on average. There- 
fore, by adding up the masses of particles with v$ < 
(and multiplying by 2) we can get an estimate of the to- 
tal mass in a spherical component. This will certainly 
underestimate /bulge in the case of rotating bulges but 
it is a reasonable approximation for the purposes of this 
work. 

Figure [3] shows that the gas component is more ro- 
tationally supported than the stellar component at all 
times, since the infalling gas can cool down and form a 
disk even if only a fraction of the initial angular momen- 
tum is retained. In the absence of strong perturbations, 
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Figure 3. Evolution of galaxy properties evaluated within the inner kpc for the stellar (red line) and gas (blue line) components of galaxies 
g54 (left) and g222 (right), two of the most massive galaxies in the sample. From top to bottom: Total mass, fraction of kinetic energy in 
ordered rotation (calculated with respect to the total angular momentum of the galaxy), bulge mass fraction (mass fraction in a spherical 
component calculated from kinematic decomposition), and specific star formation rate. 



the stellar K rot tends to increase as the stars form from 
the gas disk (e.g., galaxy g54 from z = 6 — > 3). Galaxy 
interactions and mergers result in gas inflows towards 
the center of galaxies and leave their imprint in K rot : an 
increase in the total gas mass within Rq usually corre- 
lates with a decrease in K rot (e.g., galaxy g222 at times 
~ 1500 and 1900 Myr). As expected, /bulge anticorrelates 
with fCrot; since both quantities are calculated based on 
particle kinematics. We see that galaxy interactions re- 
sult in a temporary increase of the total gas mass as well 
as the mass fraction in a spherical component. Corre- 
spondingly, the evolution of the specific SFR within the 
inner kiloparsec appears to be modulated by galaxy in- 
teractions and mergers. These events generally trigger 
temporary increases in the specific SFR that correlate 
with changes in the gas morphology and kinematics. The 
morphology of the stellar component can be affected by 
this enhanced star formation activity, since stars forming 
from low angular momentum gas after a merger can re- 
sult in the increase of the stellar bulge component (e.g., 
galaxy g222 at t ~ 1500 Myr). 

In summary, global properties of simulated galax- 
ies and their evolution on cosmological time scales are 
dominated by smooth accretion and wi nd-regulated star 
formation ([Angles- Alcazar et al.l I2013D . Despite this, 
galaxy interactions and mergers can have a significant 
impact on the morphology and star formation properties 
at the inner regions of galaxies. As we show in the next 
section, these events will have an effect on the inferred 
central black hole accretion rates and the evolution of 
AGN activity. 



4. CENTRAL ACCRETION RATES 



We infer the accretion rates onto black holes located at 
the center of each simulated galaxy based on two different 
accretion prescriptions. The first is the spherical Bondi 
accretion, in which the mass accretion rate is determined 
by the amount of matter captured gravitationally by the 
black hole. The second is the accretion rate driven by 
torques produced by gravitational instabilities at differ- 
ent scales in the in ner galaxy. We utiliz e the analytic 
model developed by IHopkins k, Quataertl (|2011[ ) to cal- 
culate the latter accretion rate. Since the evolution of 
the central kiloparsec of galaxies has been fully charac- 
terized, we can evaluate black hole accretion rates result- 
ing from different physical processes and investigate the 
connection of AGN activity and host galaxy properties. 
As we show below, the use of a particular black hole ac- 
cretion model in galaxy formation simulations can lead 
to rather different predictions on the galaxy-AGN con- 
nection. 

We choose the most bound particle in each galaxy to 
represent the central black hole. However, in this set of 
simulations, we do not sclf-consistently update the black 
hole mass and, therefore, we neglect the gravitational 
influence of the black hole at t he scales reso lved in the 
simulation. The Bondi radius (|Bondi|[l952D of a black 
hole of mass Mbh = 10 7 M Q , assuming a sound speed 
c s = SOkms^ 1 , is re = GA/bh/Cs ~ 50 pc, which is 
below the physical softening length in our simulations for 
the redshift range of interest (e « 75-224 pc at z = 8-2). 
Furthermore, the typical masses of central black holes 
that we infer for our simulated galaxies (see Section [5]) 
are, at most, comparable to the mass of a few tens of gas 
particles and we typically resolve the inner kiloparsec of 
galaxies with thousands of gas particles. Therefore, given 
the mass and force resolution in our simulations, we do 
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Figure 4. Inferred Bondi rates for galaxies g54 (left) and g222 (right). Top: Gas density at the location of the most bound gas particle. 
Middle: Sound speed at the location of the most bound gas particle (solid line) and relative velocity with respect to the surrounding gas 
(dotted line). Bottom: Inferred black hole accretion rates using the Bondi-Hoyle-Littleton parametrization (Equation [l] with a = 100) and 
based on the gas properties at the location of the most bound gas particle, where we have used a constant black hole mass Mbh = 10 Mq at 
all times and normalized to the peak accretion rate. 




Figure 5. Inferred gravitational torque rates for galaxies g54 (left) and g222 (right). From top to bottom: (i) total (stellar and gas) 
disk mass evaluated within Ro = lkpc, (ii) total disk mass fraction within Ro, (Hi) ratio of gas mass to total (stellar and gas) disk mass 
evaluated at Ro (solid line); provided that / ga B > fo (dashed line) inflow rates are not limited by gas supply, and (iv) inferred black 
hole accretion rates using the analytic model of Hopkins & Quatacrt (2011) (Equation where we have used a constant black hole mass 
Afj3H = 10 5 Mq at all times and normalized to the peak accretion rate. 
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not expect the lack of gravitational force from the central 
black hole to have a significant effect on accretion rate 
estimates. It can be argued, however, that a more critical 
issue is the lack of black hole feedback in our simulations. 
We will address the issue of feedback in Section [5l 

4.1. Bondi-Hoyle- Littleton Parametrization 

The Bondi model dHoyle fe Lvttletonl 119391 : 
iBondi fc Hovlelll944l ; lBondilll952t) is the most widely used 
prescription for b lack hole growth in galaxy formation 
simul ations (e.g.. iSpringel et all 120051 : iDi Matteo et al.l 
120081 ). For a black hole of mass .MbHj moving at velocity 
v relative to a uniform distribution of gas with density 
p and sound speed c s , the Bondi rate is given by 



M, 



Bondi 



4nG 2 M| H p 

( c 2 + ^2)3/2 ' 



(1) 



where a is a dimensionlcss parameter which is usu- 
ally added to boost accretion rates and partially com- 
pensate for the relatively high mean gas temperatures 
resulting from the multi-phase sub-grid model of star 
formation and/or the lack of the spat ial resolution re- 
quired to resolve the Bond i radius (jBooth fc Schayd 
120091: Uohansson et al.l I2009D . The choice of this pa- 
rameter, together with the initial black hole mass, can 
have a significant effect on the early growth of black 
holes. Here we use a constant v alue a = 100, sim i- 
lar to many previous studies (e.g.. ISpringel et al~ll2005j ). 
For compari son, we also explore the functional form in- 
troduced by iBooth fc Schavel (|2009[ ) , where a oc p 2 for 
gas at densities above the threshold for star formation 
(n > 0.13 cm -3 in our simulations) and a — 1 for lower 
density, single-phase gas. We note that further modi- 
fications to the Bondi parameterization have been pro- 
posed for the case of efficient cooling and significant con- 
tribution of the surrounding halo to the total gravita- 
tional potential (|Hobbs et al.l 120121 ). Making the reason- 
able assumption that the black hole is located at the 
center of the potential well, we can get an estimate of 
the Bondi rate based on the properties (gas density and 
sound speed) of the most bound gas particle at the center 
of each simulated galaxy. 

Figure 0] shows the evolution of the density and sound 
speed at the location of the black hole in galaxies g54 and 
g222. Beginning at z = 8, the central density of galaxy 
g54 increases by over two orders of magnitude up to ~ 
250 cm -3 at the end of the last merger (t ~ 2000 Myr), 
and then decreases rapidly down to ~ 30 cm -3 at z = 
2. Prior to the last merger, the accretion of low mass 
gas-rich satellites result in temporary increases of the 
central density that correlate with morphological changes 
and increases in specific SFRs as seen in Figure [3] The 
effects of mergers in the central density of galaxy g222 
are more evident, with significant density peaks (up to ~ 
400 cm -3 ) clearly correlating with changes in the specific 
SFR within the inner kiloparsec. 

We evaluate the relative velocity of the surrounding 
gas (v) as a SPH-kernel weighted averaged with respect 
to the most bound gas particle, resulting in significantly 
lower values compared to the typical sound speed (Fig- 
ure 2]). For both galaxies, the evolution of the central 
sound speed resembles that of the density but with less 
than a factor of 3 variation during the simulation. We 



note that the large values of sound speed shown here are 
due to the effective equation of state resulting from the 
sub-grid prescription of star formation. This leads to 
a significant suppression of Bondi rates that is partially 
compensated b y the addition of the boost factor a in 
Equation [T1 ( see IPelupessv et al.ll2007l for an alternative 
approach) . 

The Bondi rates inferred for black holes located at the 
centers of galaxies g54 and g222 are shown in Figure |4j 
Since the goal here is to identify the physical drivers of 
black hole accretion, we simply evaluate Equation [T] for 
a nominal, constant black hole mass Mbh = 10 5 M at 
all times and show normalized Bondi rates relative to 
the peak value for each galaxy. This simplification al- 
lows us to avoid making assumptions about the actual 
growth of black holes at this stage and to focus only on 
how the evolution of the host galaxy relates to relative 
changes in black hole accretion according to the Bondi 
model. Figure H] shows that the evolution of the normal- 
ized Bondi rates resembles that of the central density 
and it is therefore imprinted with the effects of galaxy 
mergers. Since mergers invariably cause an increase in 
the central density of galaxies, the Bondi model predicts 
a direct connection between AGN activity and galaxy 
merg ers, as has been found i n many previous sim ulations 
(e.g.. IDi Matteo et al.ll2005l: ISpringel et al.ll2005[) . More- 
over, since star formation is also proportional to the gas 
density (the sub- grid model is tu ned to match the ob- 
served relation of iKennicuttll 1998b . the Bondi model pro- 
vides support for a merger-drive n scenario for the ori gin 
of starbursts and quasar phases (jHopkins et al.|[2006f) . 

4.2. Accretion due to Gravitational Torques 

In any realistic flow with a non-zero net angular mo- 
mentum, the gas is expected to settle down onto a rota- 
tionally supported disk. The rate at which the gas is ac- 
creted depends on the rate at which its angular momen- 
tum is removed. Here, we evaluate black hole accretion 
rates b ased on the analytic model of lHopkins fc Quataertl 
(|2011[ ) that accounts for the angular momentum trans- 
port in the galactic disk. 

The gravitational torque model predicts inflow rates 
at sub-parsec scales as a function of physical quantities 
evaluated at scales that can be resolved in the simula- 
tion. Based on this prescription, the black hole accre- 
tion rate (i) increases linearly with the total (stellar and 
gas) disk mass, Mdisk, («*) depends strongly on the total 
(stellar and gas) disk mass fraction, /disk, (Hi) increases 
with the ratio of gas mass to total disk mass, / gas , and 
iv) depends very weakly o n the black hole mass, Mbh 
Hopkins fc Quataertll2011f ): 

,-, ,5/2 ( M BH A 1/6 / M disk (i? )V 

M Torquc « a T / disk x I J {-j^ 



Ro 
100 pc 



-3/2 



fa 



gas 



M yr -1 , 



(2) 



where 

fo « 0.31 / d 2 isk (M disk (i?o)/1O 9 M )- 1 / 3 , (3) 

/ gas (i?o) = M gas (i? )/Mdi sk (i?o), (4) 
and ax is a normalization factor of order ~ 1-10 that 
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parameterizes the dependence of inflow rates on star for- 
mation at scales not resol ved. Here we use cut = 5 
(|Hopkins fc Quataertl l201lh and evaluate all quantities 
within a radius Rq = 1 kpc that is well resolved in our 
simulations. The total disk mass (A/disk) is calculated 
from kinematic decomposition (similar to /bulge in Sec- 
tion 13. 2p and the total disk mass fraction is calculated 
as /di sk = M disk /(M gas (i? )+A/ s tar(-Ro)), with M gas (i? ) 
and M s tar(-Ro) the total gas and stellar masses within Rq 
respectively. 

Figure [5] shows the evolution of Mdi sk , /disk, and / gas 
for galaxies g54 and g222 as well as the resulting black 
hole accretion rates according to the gravitational torque 
model (evaluated for a constant black hole mass Mbh = 
10 5 M Q , as before). In our simulations, relative changes 
in the gravitational torque rate primarily follow from the 
evolution of Md; s k and are, therefore, very sensitive to 
the morphology of the inner region of the galaxy. The 
ratio of the gas mass to the total disk mass, / gas , tends 
to decrease with time as gas is transformed into stars 
and the disk mass increases, except during galaxy inter- 
actions that cause a temporary increase in gas mass and 
bulge fraction within Rq (e.g., galaxy g222 at ~ 1500, 
1900 Myr). Provided that / gas > /o, as it is the case in 
our simulations down to z = 2, the gravitational torque 
rate is fairly insensitive to the gas fraction, since densi- 
ties at small scales are set b y an equilibrium between g as 
inflows and star formation (|Hopkins fc Quataertll2011l ). 

The total disk mass within the inner kiloparsec of 
galaxy g54 increases by a factor of ~ 5 from z = 6 to 
z = 3. Most of the accreted gas settles onto a rota- 
tionally supported disk from which stars form, resulting 
in a significant decrease in the stellar mass fraction in 
a spherical component (/bulge; Figure [3]). Interestingly, 
galaxy mergers between z sw 3-6 seem to increase the disk 
mass and therefore the black hole accretion rate, except 
for a temporary decrease in disk mass at the end of the 
last merger (t ~ 2000 Myr) that results in a reduction of 
black hole accretion that lasts ~ 100 Myr, until the disk 
fraction rises again. 

So far, we inferred accretion rates for a constant black 
hole mass at all times and it is, therefore, not possible 
to provide a direct comparison between the Bondi and 
gravitational torque models. Despite this, the relative 
changes in black hole accretion as predicted by the two 
models seem roughly consistent for galaxy g54, with sig- 
nificant increases in black hole accretion rates correlat- 
ing with galaxy mergers. However, this scenario changes 
when looking at the evolution of galaxy g222: mergers at 
times t ~ 1500 and 1900 Myr cause a significant increase 
in Bondi rates (due to the increase in central density; Fig- 
ure 3]) but not in the gravitational torque rates. In this 
case, galaxy mergers inhibit the formation of a massive 
disk component within the inner kiloparsec. The bulge 
gas fraction increases significantly for a period of ~ 100- 
200 Myr after each merger, from which stars quickly form 
without contributing to the disk component (Figure [3]). 
It is not until the end of the second merger that the disk 
mass increases for both the gas and stellar components 
and the black hole accretion rate increases according to 
the gravitational torque model. Interestingly, the last 
two mergers at times t ~ 2700 and 3000 Myr favor the 
rapid formation of a massive disk component with the 



consequent increase in black hole accretion rates. 

In summary, galaxy mergers always cause an increase 
in the Bondi rate due to the increase in density at the 
center of galaxies, but the effects of mergers on the grav- 
itational torque rates depend on whether they cause an 
increase or a decrease in the disk mass fraction. Large- 
scale gravitational torques produced by mergers have 
the potential to remove angular momentum and drive 
gas inflows towards the centers of galaxies, as we have 
seen for two of our simulated galaxies. However, if a 
merger results in the overall reduction of the disk com- 
ponent, further gravitational instabilities in the disk are 
suppressed, with the consequent reduction in gas in- 
flow rates. The gravitational torque model yields, there- 
fore, a less direct correspondence between major merger 
events and rapid phases of black hole growth, which 
may help explain observational studies showing that ac- 
tive galaxies do not preferentially show merger signa- 
tures a t z ~ 2 (iKoceyski et al.l l2012t iMullanev et al.l 
l2012bt [Schawinski et al. 20121) and at lower redshifts 
(iGabor et alj 120091: iCisternas et alj|2011at iBoehm et al.l 
120121 ). 

5. CONSTRAINTS ON BLACK HOLE GROWTH 

Accretion rates as predicted by the Bondi and gravita- 
tional torque models depend on the effects of galaxy in- 
teractions and mergers on the evolution of central galax- 
ies. More importantly, these two models have a very 
different dependence on black hole mass. Thus, in or- 
der to provide a more meaningful comparison between 
these models and get a better intuition on the physical 
processes governing cosmological black hole growth they 
should be compared for a black hole that is growing with 
time in some consistent way. 

If we assume that the observed z = black hole-galaxy 
correlations are universal and do not depend strongly on 
rcdshift, then we can get a reasonable estimate of black 
hole masses in terms of the properties of each simulated 
galaxy. Recent observations of active galaxies indicate 
a possible evolution of the MRH- Mhnip-R (and perhaps 
-ATrh-ct) relation w ith redshift (e.g.. iDecarli et aDl2010b 
iMerloni et all 120101 ) . Therefore, black hole masses de- 
rived from the black hole-galaxy correlations should be 
taken just as a convenient parametrization as a function 
of time. In this way, we can infer the absolute instanta- 
neous accretion rates predicted by the Bondi and gravi- 
tational torque models and evaluate the implications of 
different physical processes in the evolution of black holes 
over cosmic time. 

Figure O shows the mass of the central black hole as a 
funct ion of time calculated according to the M^-a rela- 
tion dTremaine ct al. 20021) and the MBH-Afbuige relation 
(|Haring fc Rixll2004f ) for all simulated galaxies. Here, we 
calculate the effective radius of each galaxy (i? e ff) as the 
two-dimensional projected radius enclosing one half of 
the total stellar mass. The galaxy velocity dispersion (a) 
is, then, evaluated as the one-dimensional stellar velocity 
dispersion within i? c ff • Both i? g and a are averaged over 
100 random lines of sight. The bulge mass is calculated 
as the stellar mass in a spherical component within i? ff 
(from kinematic decomposition; see Section [3. 2 p and also 
taken as the total stellar mass within R c g for simplicity. 
Black hole masses shown in Figure [6] have been averaged 
over time intervals of ~ 200 Myr. 
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Figure 6. Inferred masses of central black holes as a function of time for all simulated galaxies assuming that the observed z = black hole- 
galax y correlations hold at all times. For each galaxy, black hole masses are calculated according to (i) the M^n—cr relation ( Trcmai ne et al.l 
2002) for the stellar velocity dispersion within the effective radius (dotted line) and (ii) the MeH-Afbuigc relation (Harms & Rix 2004) for 
the total stellar mass (solid line) and the bulge mass (dashed line) within the effective radius. For each galaxy, black hole masses have been 
averaged over time scales of ~ 200 Myr. 




Figure 7. Inferred accretion rates as a function of time for a black hole located at the center of each simulated galaxy. Black hole masses 
are taken from the AfeH-A^bulge relation for the total stellar mass within the effective radius at all times (FigureUl solid line) and accretion 
rates are calculated based on the properties of the host galaxy over time, according to (i) t he Bo ndi-Hoyle-Littleton parametrization with 
a = 100 (solid b lue line) and for the d ensity -dependent a introduced by Booth & Schaye (2009) (dashed blue line), (ii) the gravitational 
torque model of Hopkins & Quatacrt (2011) (red line), and (Hi) the Eddington rate (black line). Dashed green lines correspond to the 
actual accretion rates required for black holes to grow according to the MeH-Afbulge relation for each galaxy at all times (i.e., the time 
derivative of black hole masses shown in Figure 16). 
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Black hole masses derived from the Mbh-c relation are 
particularly noisy even after smoothing over long time 
scales due to the finite resolution and the calculation of a 
from individual particle motions. We note that the stel- 
lar bulge masses inferred from kinematic decomposition 
usually represent lower limits to the "true" bulge mass 
and it is therefore not surprising that black hole masses 
decrease with time according to the MBH _ -^buige rela- 
tion during certain periods of galaxy evolution. Despite 
this and given the inherently complex evolution of sim- 
ulated galaxies, it is encouraging that the Mbh _ c and 
MBH _ -^buige relations predict black hole masses that are 
largely consistent with each other. 

We can now compare the Bondi and gravitational 
torque models over cosmic time by re-evaluating Equa- 
tions [T] and [5] for a black hole mass growing according to 
the black hole-galaxy correlations. Figure [7] shows the 
evolution of the Bondi and gravitational torque rates, as 
well as the corresponding Eddington rates, for all sim- 
ulated galaxies. For simplicity, black hole masses have 
been taken from the Mbh - -Mbuige relation for the total 
stellar mass within R c e (solid lines in Figure [6|). We 
note that black hole masses are fixed by the stellar mass 
of the parent galaxy over time and therefore the inferred 
Bondi, gravitational torque, and Eddington rates shown 
here represent instantaneous accretion rates at a given 
time and do not reflect the actual growth of black holes 
according to each model. 

Some general trends can be identified in Figure [7l At 
very early times, black holes are small (with masses in the 
range ~ 10 3 -10 5 M according to the MsH-Afbuige rela- 
tion of simulated galaxies) and this results in a signif- 
icant suppression of Bondi rates due to the strong de- 
pendence on black hole mass (oc Mg H ). The Eddington 
rate increases linearly with black hole mass and, there- 
fore, the inferred instantaneous Bondi rate may eventu- 
ally become significantly higher than the corresponding 
Eddington limit for sufficiently massive black holes. In 
contrast, the gravitational torque model has a very weak 

dependence on black hole mass (oc and predicts 

gas inflows from galactic scales to sub-parsec scales that 
can exceed the Eddington limit by an order of magnitude 
at very early epochs. The gravitational torque rate also 
increases with black hole mass but is significantly more 
sensitive to the evolution of galaxy properties. Accord- 
ing to this model, black holes would experience phases of 
both super-Eddington and sub-Eddington feeding from 
z = 8 to z = 2. 

The dependence of accretion models on black hole 
mass can have profound consequences on the inferred 
black hole evolution. If black hole growth is limited 
by the gravitational capture of gas as in Bondi, black 
holes will accrete gas at very low rates early on and may 
never reach the conditions fo r critical Eddington growth 
(e.g., iPelupessv et al.l [20071. The rapid early growth 
of supermassive black holes imp l ied by observations o f 
z > 6 quasars (jFan et al.l [2005 iMortlock etldl 1201 1[ ). 
then, requires the formation of massive black hole seeds 
(~ 1 5 M ( 7 1 ) for which early accretion is less suppressed 
(e.g. lDi Matteo et al.ll20"o1l . Once black holes are mas- 
sive enough, the Bondi model results in accretion rates 
well above the Eddington rate and an additional mecha- 
nism regulating black hole growth is required. 



The black hole mass at which the Bondi rate transi- 
tions from being sub-Eddington to super-Eddington de- 
pends on mod el parameters and technical implementa - 
tions (see, e.g.. lBooth fc Schavd[200l iChoi et al.ll2012| ). 
Wc illustrate this in Figure [7| by showing how the density- 
dependent boost factor a oc p 2 results in accretion rates 
that may differ by several orders of magnitude compared 
to the constant a = 100 model, shifting the transition 
from sub- to super-Eddington to earlier times and, there- 
fore, to lower black hole masses. Despite this, the strong 
dependence of Bondi rates on the black hole mass has 
clear consequences at the low mass and high mass limits 
regardless of the particular value of a or other implemen- 
tation details. 

If black hole growth is limited primarily by the trans- 
port of angular momentum in the galactic disk, early 
growth can be rapid, since the inflow rates predicted by 
the gravitational torque model can be well above the Ed- 
dington limit even for small initial black holes. Gas in- 
flows are driven by global gravitational instabilities in 
the disk and therefore do not depend strongly on black 
hole mass. Thus, the gravitational torque model may 
ease constraints on models of black hole seed formation 
(|Volonterill20T0h . 

Figure [7] also shows the actual accretion rates required 
for black holes to grow according to the MBH-A^buige re- 
lation for each galaxy at all times (M sc i). Here we eval- 
uate M sc \ numerically as the time derivative of the black 
hole mass calculated from the MBH _ -Wbuige relation for 
each galaxy as a function of time (Figure [6]) . We find 
that the inferred Bondi rates are significantly lower than 
Msd at very early times and become significantly higher 
than M sc \ towards z = 2 for all simulated galaxies. This 
is a consequence of the strong dependence on black hole 
mass, as discussed above. In contrast, the gravitational 
torque rate seems to follow the same dependence as M sc \, 
but it is about an order of magnitude higher at all times. 
If black hole growth is limited primarily by the trans- 
port of angular momentum by gravitational instabilities, 
it is, therefore, not expected that all of the mass that is 
fed into the accretion flow at < 0.01 pc scales finds its 
way down to the central black hole. In fact, a number 
of theoretical and observational studies have shown that 
only a small fraction of the mass inflow is retained in the 
accre tion flow, with the rest lost to winds and outflows 
(e.g. JDi Matteo et al.lfeOOOt King et al.ll2012t lYuan et al.l 
12012 ). 

We compare in Figure [8] the inferred Bondi rates (with 
a = 100) and gravitational torque rates to M S c\ for all 
galaxies at all times. Remarkably, a simple power law 
can provide a reasonable fit for all simulated galaxies. 
Specifically, wc find 

Meondi = (6.07 ± O.Oe)^ 7 ^ ^ (5) 

M Torquc = (1.26 ± 0.03)^ 004±0 - 01 °) (6) 

We note that, in most cases, large deviations from 
the best fit are due to the calculation of black hole 
masses (and therefore M sc \) directly from the evolution 
of galaxy masses. The large power index (~ 2.7) for 
the Bondi rate confirms the trends found in Figure [7] 
for all simulated galaxies: Bondi rates evolve from be- 
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Figure 8. Bondi rates with a = 100 (top) and gravitational torque rates (bottom) as a function of the accretion rate required for black 
holes to grow according to the Mg[[- Mt, u i ge relation for each galaxy at all times (M sc \). Each color corresponds to a different galaxy as 
in Figure [2] Solid lines show the best power law fit for each model and dashed-dotted lines correspond to a linear relation with constant 
slope given by median(M^ on< i L [ / M sc \) (top) and median(M^ OVCL ^ e / M sc i) (bottom) for all galaxies at all times. 



ing significantly lower than M sc i at small black hole 
masses to being significantly higher at large black hole 
masses. That is, black holes growing at the Bondi 
rate would either lie well below the -MBH~-^buige rela- 
tion at all times or would become overly massive by 
z = 2. If the Bondi rate provides a good estimate of 
the true accretion rate, the large power index in Equa- 
tion[5]suggests that (i) the initial black hole mass and the 
boost factor a need to be set to allow for early growth 
and (it) some form of feedback that becomes more effi- 
cient at higher accretion rates is required. Indeed, suit- 
able choices of model parameters and the coupling of a 
fraction of the accretion luminosity to the surrounding 
gas in the form of thermal energy have been success- 
ful in rep roducing the observed black hole-galaxy cor- 
relations (iDi Matteo et all 120051: iRobertson et al.l 12006b 
IHopkins et alJl2007t IBooth fc Schavell2009D . 

The gravitational torque rates are, on average, di- 
rectly proportional to the accretion rates required by 
the MBH _ -Mbuigc relation, with a power index very 
close to unity (Equation [6]). Indeed, a linear relation 
-^Torque = e^ 1 Msci with constant slope given by e^ 1 = 
median(MT 0vque / M sc \) for all galaxies at all times, pro- 
vides a very reasonable fit to the numerical results (Fig- 
ure [5]). Remarkably, a simple mass retention rate in the 
accretion flow, e m rj 5%, seems sufficient to bring the 
gravitational torque rates down to M sc \ over a range of 
more than three orders of magnitude in accretion rates, 
for all simulated galaxies, and at all times. Therefore, 
the additional mechanism required to regulate black hole 
growth according to the MBR-M^ u i ge relation has to be 
roughly equally efficient regardless of the accretion rate. 



Since the radiative luminosity is proportional to the black 
hole accretion rate, the gravitational torque model seems 
to disfavor AGN feedback acting at galactic scales as the 
primary mechanism regulating black hole growth. 

6. TORQUE- LIMITED GROWTH 

The transport of angular momentum in the galactic 
disk is a required process for black hole growth regardless 
of the presence and effects of AGN feedback. Therefore, 
we can explore an alternative scenario in which winds 
and outflows from the inner accretion disk do not affect 
the accretion flow and simply result in the reduction of 
black hole accretion rates by a constant mass retention 
rate e m . The gravitational torque model describes gas 
inflows from galactic scales down to < 0.1 pc scales and 
therefore it can be interpreted as the mass feeding rate 
into the black hole accretion disk. Some of this mass may 
be los t by radiatively dr i ven outflows and other processes 
fe-g JProga et al.l 120001: lOhsuga et all 120051: lYuan et al.1 
|2012|) . while a small fraction reaches the central black 
hole. In this scenario, the central black hole grows on 
average at a fraction e m of the large-scale gas inflow rate: 



dM-BH./dt = e m Mrorque(t) 



(7) 



Figure [9] shows the evolution of black hole masses pre- 
dicted by the gravitational torque model normalized by 
the mass retention rare e m 0.05 (M^ quc ), as a func- 
tion of the black hole mass obtained from the Mbh - 
-Wbuige relation for each galaxy at each time step (M B jj). 
Here, we are simply integrating Equation [7] for an ini- 
tial black hole mass which is consistent with the A/bh~ 
-Mbuige relation for each galaxy. We note that MrorqucW 
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is self-consistently calculated based on the morpholog- 
ical properties of each galaxy over time (Equation [3]) 
and updated with the appropriate black hole mass at 
each time step, as given by Equation Remarkably, 
we find that black holes evolve approximately along the 
-^bh qU ° = ^bh nne ano - therefore, remaining consistent 
with the A/bh ^buigc relation. 

In Section 17.21 we show that e m sets the normaliza- 
tion of the MBH _ -Mbulge relation. We note, however, 
that the gravitational torque model (Equation [2]) con- 
tains a normalization constant of order ctt ~ 1-10 that 
parametrizes the ef fects of nuclear star form ation on 
black hole accretion ([Hopkins fc Quataertl[20TT[ ). There- 
fore, the normalization of the Mbh - A/bulge relation is 
actually determined by the product c*t e m . Given that 
plausible values of these two constants differ by about 
two orders of magnitude, we choose to interpret them 
as parametrizing different physical processes: q?t con- 
trols how much gas is converted i nto stars versus feed- 
ing t he black hole accretion disk ([Hopkins fe Quataertl 
1201 lh and e m controls what fraction of the mass feeding 
the accretion disk is finally accreted by the black hole. 
We implicitly assume that the remaining mass, a fraction 
(1 — e m ) of the disk feeding rate, it is expelled in the polar 
direction without significantly halting further accretion. 

We can now look back at the accretion histories of 
black holes growing according to the gravitational torque 
model with the appropriate normalization (e m w 0.05). 
Figure \W\ shows the evolution of the inferred Eddington 
ratios of black holes located at the center of all simulated 
galaxies. Despite the large scatter between galaxies and 
the high temporal variability, there is a clear trend for 
Eddington ratios to decrease from ~ 0.1-1 at z = 8 down 
to ~ 0.01-0.1 at z = 2. This is in agreement with recent 
observations of z > 4 quasars showing that higher red- 
shift black holes have lower masses but at are accreting at 
higher Eddington ratios compared to lower redshift b lack 
holes (|De Rosa et al.ll2011t [Trakhtenbrot et al.ll2011ft . 

Intriguingly, the inferred Eddington ratios evolve in a 
manner very similar to the evolution of specific SFRs 
(Figure [2]): both quantities seem to decrease by about 
a factor of 10 from z = 8 — > 2, on average. In fact, 
Figure [10] (lower panel) shows that the ratio of the black 
hole accretion rate to the SFR within the effective radius 
of galaxies is, on average, very close to the ratio of black 
hole mass to stellar mass as inferred from the local A/bh- 
Afbuige relation: 

( Trll ) ~ 0.00154 ± 0.0008, (8) 
\ or K I 

where the average (and standard deviation) has been cal- 
culated for all galaxies and all times. There is, of course, 
large scatter in A/bh/SFR between different galaxies and 
evolution times, but this ratio is, on average, remarkably 
flat in the redshift range z = 2-8: the bulk of black hole 
and galaxy growth occurs in tandem, at rates governed 
by cosmological infall and transport of angular momen- 
tum in the galaxy disk, and with no feedback-mediated 
coupling between black holes and parent galaxies. 

The torque-limited model predicts, therefore, a natu- 
ral connection between AGN activity and SFR of galaxies 
on cosmological time scales, in agreement with observa- 
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Figure 9. Black hole mass calculated self-consistently from the 
gravitational torque rate corrected by the mass retention rate e m = 
0.05 (by integration of Equation[jJ as a function of black hole mass 
derived from the JWbh— -^bulge relation for each galaxy at all times. 
Each color corresponds to a different galaxy as in Figure [2] 
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Figure 10. Top: Eddington ratio (defined here as Men/A^Edd) 
for black holes growing self-consistently at the gravitational torque 
rate corrected by the mass retention rate e m = 0.05 (by integra- 
tion of Equation [jj for each galaxy as a function of time. Bottom: 
Ratio of the black hole accretion rate to the total SFR within the 
effective radius. Black solid and dashed lines show a running aver- 
age of Mbh/SFR over all galaxies and all times and the standard 
deviation respectively. Initial black hole masses are taken to be 
consistent with the 2 = M^n-M hu i se relation (Figure [BJ. Each 
color corresponds to a different galaxy as in Figure [3] 
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tions showing that the volume-averaged ratio of black 
hole accretion rate to SFR agrees with the ratio of black 
hole mass to stellar mass as inferred from the local Mbh - 
Mbuige relation (jHeckman et al.ll2004l : IZheng et all2009h . 
Furthermore, these results are consistent with recent ob- 
servational evidence for a constant black hole accretion to 
star formation ratio independent of redshif t since z ~ 2 
(jRaffertv et al.llMlUMullanev et al.ll2012al) . Since black 
holes and galaxies evolve approximately along the scal- 
ing relations, the torque-limited model provides support 
for a roughly constant black hole mass to stellar mass 
ratio inde pendent of redshift. consistent with the con - 
clusi ons oflJahnke et all (120091) , iCisternas et~aTl (|2011bl ) , 
and iMullanev et all ()2012aF at redshifts z < 2. The 
AGN-SFR connection is, however, not direct in a galaxy- 
by-galaxy basis and at all times (as evident from the scat- 
ter in Figure [lOj lower panel), which may help explain the 
large scatter in observed AGN to star formati on ratios 
in individual systems (e.g.. iRaffertv et aT1l2011l) . 

Torque-limited growth is, therefore, a natural alter- 
native to feedback-regulated models and it is indeed 
consistent with our "feedback- free" starting point, since 
black hole accretion rates have been inferred based on 
galaxy properties from simulations with no AGN feed- 
back model (Section 3]). AGN feedback may have a sig- 
nificant impact on the host galaxy but it is not required 
for regulating black hole growth. Angular momentum 
transport (in the galaxy as well as the accretion disk) to- 
gether with competition with star formation seem to be 
sufficient for black holes to grow according to the scal- 
ing relations from early times down to z = 2. This is 
explored in more detail in the next section. 

7. BLACK HOLE-GALAXY CORRELATIONS 

Figure [TT] shows how galaxies and black holes evolve 
in the MBH _ -Wbuige and Mbh-u planes according to the 
torque- limited model (Equation [7]). As we have seen 
in the previous section, a simple normalization con- 
stant e m ~ 0.05 applied to the gravitational torque 
model brings blac k holes and galaxies c lose to the A/bh~ 
Mhuigc relation of iHaring fc Rbd (|2004| ). Calculation of 
the galaxy stellar velocity dispersion from particle mo- 
tions results in higher scatter but still shows that black 
holes growing according to Equ ation[TJare roughly consis- 
tent with the Mbh-ct relation of lMcConnell et al.1 (|2011l ). 

In Section[5]we have inferred "instantaneous" accretion 
rates according to the Bondi and gravitational torque 
models by assuming that the observed MBH^-Mbuigc re ~ 
lation holds at all times. This was a convenient way 
of parametrizing black hole mass as a function of time 
for each galaxy, but there is a priori no reason to think 
that galaxies and black holes should evolve together at 
all times. The gravitational torque rates inferred in this 
way were often of the order of the Eddington limit (Fig- 
ure [7]) and therefore we might be tempted to think that 
black holes growing at a fixed fraction of the Edding- 
ton limit could also grow in a manner consistent with 
the black hole-galaxy correlations. However, once we 
release the assumption of black hole-galaxy coevolution 
(so that black hole masses are no longer given by the 
MBH-Afbuigo relation) and calculate black hole growth 
sclf-consistcntly over time, we see that a simple nor- 
malization constant cannot compensate for the exponen- 
tial Eddington growth. The fact that the torque-limited 



model naturally leads to black holes growing on average 
according to the black hole-galaxy correlations is a non- 
trivial consequence of the weak dependence on black hole 
mass and the strong dependence on galaxy scale proper- 
ties. 

In the remainder of this section, we show how the in- 
ferred black hole-galaxy correlations depend on model 
parameters and implementation details, and discuss ad- 
ditional implications of the torque-limited model. 

7.1. Initial Black Hole Mass 

We have shown that provided the initial black hole 
mass lies approximately in the -MBH - -^buige relation, the 
gravitational torque model scaled by a simple constant 
factor predicts that black holes and galaxies evolve to- 
gether along the scaling relations. However, black holes 
do not necessarily know what is the correct initial mass 
for a given galaxy. In Figure [T2l we evaluate how changes 
in the initial conditions affect black hole growth and the 
inferred MeH _ Afbuige and Mbh _ c relations. First, we 
grow black holes with initial masses either a factor of 10 
above or below the scaling relations for each galaxy. We 
find that, regardless of the initial mass, black holes tend 
to evolve onto the observed black hole-galaxy correla- 
tions. In our simulations, accretion rates are governed 
by the disk mass, with a very weak dependence on black 
hole mass, and therefore black holes with different initial 
masses grow at comparatively similar rates for a given 
galaxy. Since most of the black hole mass comes from 
accretion, which depends mostly on the evolution of the 
galaxy, the initial conditions are erased and the black 
hole-galaxy correlations are established. 

This is also observed when the initial black hole masses 
are uncorrelated with their parent galaxies. Figure [12] 
(middle panels) shows the -MBH-Mbuige and Mbh^ct re- 
lations when initial black hole masses are taken to be the 
same for all galaxies and ranging from 10 2 to 10 5 M Q . 
For simplicity, we show the evolution of the "average 
galaxy" that we calculate by averaging black hole mass, 
bulge mass, and velocity dispersion over all galaxies at 
each time step. Provided we let black holes and galax- 
ies evolve for a sufficient amount of time, they all con- 
verge on average towards the scaling relations. Finally, 
we show in Figure [12] (bottom panels) that this conver- 
gence towards the scaling relations is also observed when 
we take black hole masses for all galaxies at all times and 
average them within bins in either bulge mass or velocity 
dispersion regardless of the evolution time of each galaxy. 

Black hole-galaxy correlations arising from the torque- 
limited model are robust to changes in the initial mass 
of black holes and do not require specific tuning of the 
mass retention rate e m . Importantly, black holes with 
masses as low as 100 M© at z = 8 are able to grow com- 
paratively faster than higher mass black holes, opening 
up possibilities for light seed formation mechanisms, such 
as remnants of population III stars, of be ing the progeni- 
tors of today's superm assive black holes (jMadau fc Reesl 
[20??ltlVolonter1 1201(1 . 

7.2. Mass Retention Rate 

The mass retention rate (e m ) is a free parameter in 
the torque-limited model that we have introduced in or- 
der to match the observed MBH-Mb n \ ge relation (Sec- 
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Figure 11. Mbh _ A^bulge and Mbh— o" (right) relations for black holes growing self-consistently at the gravitational torque rate 

corrected by the mass retention rate e m = 0.05 (by integration of Equation for each galaxy at all times. Initial black h ole masses are 
taken to be already consistent with the scaling relations. Black solid lines show the Mg]j- -Mbui ge and Mbh-it relations of IHarine &: Rixl 



(2004) and McConncll ct al. (2011) respectively and black dashed lines indicate a 0.5 dex scatter in black hole mass for each relation. The 
galaxy velocity dispersion (<r) is calculated as the one-dimensional stellar velocity dispersion within the effective radius (-R c ff ) averaged over 
100 random lines of sight. Mt,ui go is taken as the total stellar mass within i? e ff. Each color corresponds to a different galaxy as in Figure |21 



tions [5] and [6]) . Figure [13] shows how the inferred A/bh - 
-^buigc relation changes when using a wide range of re- 
tention values, from a extreme situation in which 99.9 % 
of the inflow mass is lost in outflows, to the upper limit 
e m = 1 in which all of the infalling gas from larger scales 
is accreted by the black hole (which corresponds to the 
original normalization of lHopkins fc Q uataort l2~01lD . As 
we might expect, we find that e m only affects the normal- 
ization of the MBH _ Afbuige relation, with little effects on 
the slope except at very early times. 

A somew hat simi l ar tor que-limited scenario has been 
proposed bv lEscalal (|2006|) and tested with idealized sub- 
parsec resolut ion simulations of a nuclear galactic disk 
(|Escalall2007t) . They find that a mass retention rate 
e m ~ 15 % is required in order to match the A/bh~ 
A^buigc relation, similar to our result e m ~ 5 % given that 
both sets of simulations probe spatial and time scales 
that are substantially different. 

7.3. Stellar Feedback and Resolution Effects 

In a separate work using the same set of simulations, 
we have shown that strong stellar feedback in the form 
of galactic outflows can have a significant impact on the 
star formation, morphologi cal, and kinematic propertie s 
of high-redshift galaxies (pVngles- Alcazar et al.l 120131 ). 
Feedback from sta r formation could therefore affect black 
hole growth (e.g.. ICenl 120121 ) and possibly the inferred 
black hole-galaxy correlations. 

Figure [T4l shows the MBH _ Afbuige relation obtained for 
the same sample of galaxies but this time from simula- 
tions that do not include models for galactic winds. We 
find that the scaling relations are still reproduced by the 
torque- limited model and perhaps even with smaller scat- 
ter. The convergence towards the MBH~Afbui g e relation 
for different initial black hole masses is retained (Sec- 
tion 17. 1|) and similar results are obtained for the AIbw 
a relation. Predictions of the gravitational torque model 
are therefore robust to changes in stellar feedback. Even 
if galaxies from simulations with no winds become sig- 
nifican tly more massive and less g as rich owing to higher 
SFRs (| Angles- Alcazar et al.ll2013f ) . black holes still grow 
according to the scaling relations in the redshift range 



z = 8 — > 2. This suggests that galaxies regulate them- 
selves via an equilibrium between inflows and outflows 
and black holes grow according to the properties of their 
host galaxies, at a rate given by the gas supply from 
galactic scales. 

We have also investigated the effects of numerical res- 
olution on the scaling relations. Figure [15] shows the 
AfBH - A/buigc relation obtained for the same sample of 
galaxies from simulations carried out with a factor of 2 
lower spatial resolution and a factor of 8 lower mass res- 
olution (including galactic outflows). Since we require 
galaxies to contain at least 100 gas particles and 100 star 
particles within Rq = 1 kpc in order to resolve their cen- 
tral morphology, these lower resolution simulations can 
follow the evolution of black holes only after their parent 
galaxies are a factor of <~ 8 more massive compared to the 
high resolution simulations. Despite the reduced resolu- 
tion, the Mbh -A/buigc relation arising from the gravita- 
tional torque model shows good numerical convergence. 

We note that the values of e m obtained for simulations 
with no galactic outflows and the low resolution simula- 
tions are slightly different from our fiducial runs but they 
are all consistent with a typical mass retention rate of a 
few percent of the gas supply. 

7.4. Black Hole Mergers 

It has been shown recently that black hole-galaxy cor- 
relations may be a natural consequence of hierarchical 
structure formation, so that an initially uncorrclated 
distribution of black hole and stellar masses tends to 
converge towards the scaling relations thro ugh succes- 
sive mergings of black holes and galaxies (jPengl 120071 : 
Uahnke fc Maccioll201lh . In this scenario, AGN feedback 
and self-regulated growth are not required, but it is as- 
sumed that galaxy mergers always result in the merg- 
ing of their central black holes. Here, we have shown 
how torque-limited growth yields black hole-galaxy cor- 
relations regardless of the initial masses of black holes 
and galaxies, and also without the need for AGN feed- 
back and self-regulated growth. We have, however, left 
off the discussion the role of black hole mergers in the 
overall black hole growth. Here, we assume that most 
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Figure 12. MeH^-Mbulge {left panels) and A/bh-c (right panels) relations for black holes growing self-consistently at the gravitational 
torque rate (Equation [7] with e m = 0.05) for different initial black hole masses. Top panels: The initial black hole mass is taken to be 
either fO times higher (red) or ten times lower (blue) compared to the scaling relation for each galaxy. Middle panels: Black hole mass 
as a function of either bulge mass or velocity dispersion for the "average" galaxy (i.e., averaging Mjh, jWb u i go , and a over all galaxies 
at each time step) for initial black hole masses ranging from 10 2 to 10 5 Mq (from blue to red). Bottom panels: Average black hole mass 
within logarithmically-spaced bins in either bulge mass or velocity dispersion, using all simulated galaxies at all time steps, and for initial 
black hole masses taken from 10 2 to 10 5 Mq (fro m blue to red). Error bars show the dispersion in each bin. Black solid and dashed 
lines correspond to the observed scaling relations (Haring & Rix 2004; McConncll ct al. 201f) and a 0.5 dex scatter in black hole mass 
respectively. 



of th e black hole mass comes from accretion ([Soltanl 
Il982f ) and we simply neglect the mass contribution from 
mergers. Given that only a fracti on of mass in our 
galaxies comes fr om major mergers (jMurali et al.ll2003 
IKeres et al.ll2005|) , we do not expect a significant shift on 
the inferred black hole-galaxy correlations due to black 
hole mergers. 

7.5. Evolution of the Scaling Relations 

Recent observations of active galaxies seem to indicate 
an evolution of the MeH-^buigo (and perhaps A/bh - 
a) relation with redshift, with black holes being more 
massive for a given galaxy mass (or velocity disper- 



sion) at higher re d shifts compared t o the local relations 
(IShen et all [2008t IWolTet~al1 [20081: fPecarli et al.1 l20l(it 
iGreene et al.ll2010t iMerloni et alJl2010D. There are also 
findings of no significant evolution |jahnke et al.ll2009f ) 
and potentially even und ermassive black hole s in z ~ 2 
infrared-selected galaxies ([Shapiro et al.ll20 0lt). In either 
case, evaluations of the redshift evolution of the black 
hole scaling relatio ns may be biased by selection effects 
(|Lauer et al.ll2007[ ). 

In the context of torque-limited growth, we find no 
evidence for evolution of the scaling relations at least 
down to z = 2, since black holes and galaxies seem to 
converge towards the scaling relations regardless of their 
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Figure 13. The effects of varying the mass retention rate from e m = 0.001 to e m = 1 by factors of 10 (from blue to red) in the Mbh- 
Afbuigc relation. Left: Black hole mass as a function of bulge mass averaged over all galaxies at each time step for initial black hole 
masses consistent with the MeH^A^bulge relation. Right: Average black hole mass within logarithmically-spaced bins in bulge mass for all 
simulated galaxies, at all time steps, and taking initial black hole masses consistent with the MeH^-^bulgc relation for each galaxy. Black 
solid and dashed lines correspond to the observed M^n— M bu i gc relation (Haring & Rix 2004) and a 0.5 dex scatter in black hole mass 
respectively. 
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Figur e 14 . Effects of stellar feedback on the AfBH-A^bulge relation. Black holes grow according to the gravitational torque rate as in 
Figure 1111 but for simulations that do not include galactic winds. Left: Initial black hole masses are taken to be consistent with the 
A4BH--Mbulge relation for each galaxy. Each color corresponds to a different galaxy as in Figure [2] Right: Average black hole mass within 
logarithmically-spaced bins in bulge mass, for all simulated galaxies, at all time steps, and for initial black hole masses ranging from 
1O 2 M to 10 5 M Q (from blue to red). The mass ret ention rate e m = 0.03 has been used. Black solid and dashed lines correspond to the 
observed AfeH^-Mbulge relation ( Hiring & Rix 2004) and a 0.5 dex scatter in black hole mass respectively. 



initial masses, provided ojt and e m do not evolve with 
rcdshift and are well represented by a constant factor. 
Therefore, significant deviations at early times could in 
principle reflect the "initial conditions" for coevolution 
of black holes and galaxies (Figure [P2jl . 

We note that e m ha s been fixed h e re to match the 
M B H-M bu i ge relation of lHaring fc Rix! (|2004D for the to- 
tal stellar mass within the effective radius rather than 
the stellar bulge mass. Given that our galaxies con- 
tain a significant stellar d isk component at z = 2 
(jAngles- Alcazar et al.ll2013l ). it is plausible that there is 
significant evolution of the MBH _ -Wbuige relation with re- 
spect to bulge mass bu t not with respect to total stellar 
mass, as sugg ested bv Uahnke et all (|2009t) . If most of 
the mass of bulges and elliptical galaxies today comes 
from redistribution of the stellar bulge+disk of z = 2 
galaxies (through mergers or other processes) no further 
black hole accretion would be required to match the lo- 
cal scaling relations. We speculate that the formation of 
compact bulges would indeed decrease the gravitational 
torque rate (since it is proportional to the disk mass), 
truncating significant black hole growth at late times. 



8. SUMMARY AND CONCLUSIONS 

We have used cosmological zoom simulations of galaxy 
formation down to z = 2, together with analytic models 
of black hole accretion, to investigate the growth of mas- 
sive black holes at the centers of galaxies without making 
any prior assumptions about the effects of AGN feed- 
back. To this end, we have compared predictions from 
the spherical Bondi mode l and an accretion model dr iven 
by gravitational torques (jHopkins fe Quataert]|201ll ). 

We find that the Bondi model presents significant 
challenges due to the strong dependence of the accretion 
rates on black hole mass. Black hole growth is signifi- 
cantly suppressed for low mass black holes and becomes 
extremely rapid for sufficiently massive black holes. 
This has two main implications for simulations of galaxy 
formation: (?) the initial black hole mass (together 
with any boost factor or normalization constant) has 
to be sufficiently massive to allow for early black hole 
growth, and (n) feedback energy and/or momentum 
needs to be injected into the surrounding gas in order 
to regulate black hole growth at later times. Because 
these implications follow from the oc Mg H dependence, 
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Figure 1 5. Effects of numerical resolution on the MeH-Afbulge relation. Black holes grow according to the gravitational torque rate as 
in Figure [HJ but for simulations with a factor of 2 lower spatial resolution and a factor of 8 lower mass resolution. Lines and colors are as 
in Figure Il4l The mass retention rate e m = 0.05 has been used. 



they should apply to any (reasonable) modification of 
the original Bondi-Hoyle-Littlcton parameterization. 
Remarkably, suitable choices of model parameters 
for different implementations have been successful in 
reproducing a broad number of observ ations, including 
the b la ck hole-galaxy co rr elatio ns fe.g..lDi Matteo et all 
2001 iRobertson et all l2006t iHopkins et all 120071; 
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Booth fc Schavell2009t Uohansson et al 

In the gravitational torque model, gas inflows 
driven by global gravitational instabilities in the disk 
and, t herefore, do not depend str ongly on the black hole 
mass (Hopkins fc Quataertll20ilj) . The resulting accre- 
tion rates imply that (i) early black hole growth could 
in principle proceed at super-critical rates, since the in- 
flow rates can be well above the Eddington limit even 
for small initial black holes, and (ii) energy (and/or mo- 
mentum) from the accretion process does not need to 
couple to galaxy-scale gas in order to regulate black hole 
growth. Consistency with the black hole-galaxy correla- 
tions simply requires the assumption that only a small 
(constant) fraction of the mass inflow is retained in the 
accretion flow, with the rest lost to winds and outflows. 
We have shown that this result is insensitive to variations 
in the initial black hole mass, stellar feedback, or other 
implementation details, and the required mass retention 
rate (e m ~ 5 %) is rou ghly consistent with observational 
fe.g.JKing et all 20121 ) and theoretical expectations (e.g., 
lYuan et al.ll2012 , and references therein). However, the 
exact value of the normalization factor e m that we in- 
fer is subject to uncertainties such as variations in the 
nuclear star formation law or the exact normalization of 
the MBH-A^buigc relation. 

Mass outflows are invoked in the torque-limited model 
in order to control what fraction (e m ) of the accretion 
disk feeding rate is finally accreted by the central black 
hole (a linear effect), but not to regulate the amount of 
gas feeding the accretion disk from galactic scales. In 
this scenario, there is no coupling between inflows and 
outflows and, therefore, black hole feedback cannot shut 
down accretion. "Feedback self-regulation" does not oc- 
cur in the strict sense of the term, since the energy out- 
put by accretion onto the black hole has no direct im- 
pact on its own fuel supply. This greatly differs from the 
non-linear feedback loop required by the Bondi model, 
large-scale AGN feedback may have a significant impact 
on the host galaxy but it is not required for regulating 



black hole growth, which may instead be limited by the 
efficiency of gravitational torques in removing angular 
momentum from the gas together with competition with 
star formation. 

In addition, torque-limited growth yields a less di- 
rect correspondence between major merger events and 
enhance AGN activity, as sugges ted by recent obser- 
vations of z ~ 2 active galaxies (iKocevski et all 120121 : 
iMullanev et~aT1l2012bl : iSchawinski et al.ll2012h ~l"nstead. 
black hole and galaxy growth are governed by cosmo- 
logical infall and transport of angular momentum in the 
galactic disk, giving rise to a time-averaged connection 
between AGN activity and SFR on cosmological time 
scales. Our findings are consistent with recent obser- 
vational evidence for a constant black hole accretion to 
star formation ratio, similar to the the ratio of black hole 
mass to stellar mass as inferred from the local A/bh _ 
Mhnire relation, and independent of redshif t since z ~ 2 
(jRaffertv et allboiH IMullanev et alJl2012af >. 

The agreement of feedback-regulated accretion mod- 
els in hydrodynamic simulations to the observed black 
hole-galaxy correlations has been interpreted as (indi- 
rect) observational evidence of both AGN feedback act- 
ing at galactic scales and the self-regulated growth of 
black holes. We note, however, that the black hole- 
galaxy correlations are indeed a primary constraint for 
the s trength of such feedback models (see iDebuhr et all 
120111 for a self-regulated model that accounts for angu- 
lar momentum transport). Given that AGN feedback 
can have profound consequences on the e volution of the 
host galaxy (e.g., iDi Matteo et all 120051 ) and that the 
scaling relations arc not necessarily a consequence of 
feedback, it would be desirable to constrain AGN feed- 
back models by others means, ideally through direct 
obser vations of mass outflow rates (e.g., iFeruglio et all 



20101) and /or more physically motivated mode ls (e.g. 



zvlw and/ or more pnysicaiiy 

Ciotti fc O striker 2007 tlNovak et all l20ll 12012^ . While 



the Bondi model, appropriately tuned, does provide a 
good match to the black hole-galaxy correlations, our 
work demonstrates that the torque-limited model pro- 
vides a viable alternative that does not require self- 
regulation. 

Our findings should apply so long as the assump- 
tions implied by the gravitational torque model are met. 
In particular, it is assumed that (i) there is a signifi- 
cant stellar component that drives the gas into shocks 
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that dissipate energy and angular momentum, and (ii) 
the amplitudes of non-axisymmetric modes are large 
enough to produce such sh ocks even at ~ 10 pc scales 
([Hopkins fc Qua tacrt 20"Tlh. These scales are, of course, 
not resolved in our simulations; galaxy centers become 
quickly stellar dominated but the formation time of the 
first star particles is poorly constrained. Other pro- 
cesses such as scattering of dense gas clouds or the trans- 
port of angular momentum by supernova and/or grav- 
itational instability-driven turbulence ma y be required 
in the very high gas -reach domai n (e.g ., lEscalal f2006t 
IBournaud et~aLl l20lH; iHobbs et al] [20111) . Finally, full 
cosmological simulations with adequate resolution down 
to z = are required in order to evaluate the impli- 
cations of torque-limited growth on a statistically sig- 
nificant sample of galaxies, making connection with the 
wealth of data available in the local and low redshift 
universe over a range of galaxy masses spanning several 
orders of magnitude. 
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